GsERF1 enhances Arabidopsis thaliana aluminum tolerance through an ethylene-mediated pathway

Ethylene response factor (ERF) transcription factors constitute a subfamily of the AP2/ERF superfamily in plants and play multiple roles in plant growth and development as well as in stress responses. In this study, the GsERF1 gene from the wild soybean BW69 line (an Al-resistant Glycine soja line) was rapidly induced in response to aluminum stress. Quantitative real-time PCR (qRT–PCR) analysis showed that the GsERF1 gene maintained a constitutive expression pattern and was induced in soybean in response to aluminum stress, with increased amounts of transcripts detected in the roots. The putative GsERF1 protein, which contains an AP2 domain, was located in the nucleus and maintained transactivation activity. In addition, under AlCl3 treatment, GsERF1 overexpression increased the relative growth rate of the roots of Arabidopsis and weakened the hematoxylin staining of hairy roots. Ethylene synthesis-related genes such as ACS4, ACS5 and ACS6 were upregulated in GsERF1 transgenic lines compared with the wild type under AlCl3 treatment. Furthermore, the expression levels of stress/ABA-responsive marker genes, including ABI1, ABI2, ABI4, ABI5 and RD29B, in the GsERF1 transgenic lines were affected by AlCl3 treatment, unlike those in the wild type. Taken together, the results indicated that overexpression of GsERF1 may enhance aluminum tolerance of Arabidopsis through an ethylene-mediated pathway and/or ABA signaling pathway, the findings of which lay a foundation for breeding soybean plants tolerant to aluminum stress. Supplementary Information The online version contains supplementary material available at 10.1186/s12870-022-03625-6.


Introduction
Heavy metal toxicity, such as aluminum (Al) toxicity, is a major limiting factor for crop production worldwide [1]. When the pH of the soil is lower than 5.0, aluminum is present in an ionic form, i.e., Al 3+ , which strongly inhibits root growth and function, reducing crop yields [2]. Plant species and varieties vary widely in their ability to tolerate aluminum toxicity. Some plant species or varieties have evolved high levels of tolerance mechanisms to survive in acidic soils. Wild soybean has been growing in acidic soils in South China for a long time, and as such, there resources available that can provide tolerance, which plays an important role in improving the stress resistance of soybean [3]. Transcription factors are involved in stress responses; transcription factors from the WRKY, bZIP and NAC families have been shown to participate in the aluminum stress response and to regulate the aluminum tolerance of plants [4][5][6]. However, the involvement of the ERF transcription factor family in the aluminum stress response has not been reported. Plants employ a complex regulatory network to cope with a variety of stresses during growth and development. A variety Open Access † Lu Li and Xingang Li are equally contributed to this work. *Correspondence: hnian@scau.edu.cn; maqibin@scau.edu.cn 5 Zengcheng Teaching and Research Base, South China Agricultural University, Guangzhou, Guangdong 510642, People's Republic of China Full list of author information is available at the end of the article of plant hormones play important roles from the beginning of sensing stress signals to the response of plants to stress. Under normal circumstances, the ethylene content in plants is maintained at a low level. However, plant ethylene content changes in response to biological stress or abiotic stress. The response of ethylene production after stress stimulation is transmitted through corresponding signal transduction pathways, which can regulate downstream genes, causing a series of reactions in plant cells and an associated response to stress [7][8][9]. Previous studies have shown that when plants are subjected heavy metal toxicity, the general response involves increased production of ethylene. For example, plants increase their production of ethylene under toxic levels of cadmium (Cd), copper (Cu), iron (Fe), nickel (Ni) and zinc (Zn). Moreover, it has been found that the change in ethylene under heavy metal stress is due to the increased expression of ethylene-related biosynthesis-related genes and/or changes in the expression of ethylene-responsive genes. Regarding these changes, it has been found that the increase in ethylene during stress can have negative effects on plants. However, ethylene can alleviate the inhibition of the photosynthetic capacity of mustard under cadmium stress. These findings suggest that ethylene involves a complex two-way regulatory function under stress, which depends on its endogenous level [10][11][12][13].
ERF transcription factors (ethylene response factors) constitute a subfamily of the AP2/ERF superfamily and can be divided into three categories according to the number of AP2/ERF domains: AP2, ERF and RAV [14]. The ERF family protein members contain an AP2/ERF domain consisting of 58-60 highly conserved amino acids, which constitutes the main functional region of ERF family proteins [15]. Ethylene response factors (ERFs) not only play important roles in plant growth and development but also play very important roles in the plant response to stress [15]. Previous studies have shown that ERF family genes are involved in plant growth and development in rice, Arabidopsis and other plant species. For example, OsERF1 is constitutively expressed in different organs of rice and is upregulated by ethylene. Overexpression of OsERF1 significantly affects the growth and development of transgenic Arabidopsis by promoting the expression of the ethylene-responsive genes PDF1.2 and β-chitinase [16]. AtERF71/HRE2 can activate the expression of downstream genes by binding the motifs of GCC boxes and DRE/CRT elements, regulate the expansion of root cells and play important roles in root development [17]. Julien Pirrello found that overexpression of the Sl-ERF2 gene in transgenic tomato lines can lead to early seed germination and enhanced hypocotyl formation in dark-grown seedlings. Recently, the transcription factor ERF139 was found in poplar to regulate the expansion of xylem cells and the deposition of secondary cell walls [18].
In recent years, an increasing number of ERF family genes have been found to function in stress tolerance in plants. Under drought stress, overexpression of the rice genes OsERF71, OsERF101 and OsERF48 was shown enhance the drought resistance of rice [19][20][21][22]. Heterologous overexpression of the soybean gene GmERF3 can enhance tobacco drought resistance [23], and overexpression of AtERF019 can enhance drought resistance in Arabidopsis [24]. Overexpression of GmERF135 can enhance the salt tolerance of Arabidopsis plants under salt stress. Moreover, GmERF135 can promote the growth of transgenic hairy roots under salt stress [25]. In wheat, overexpression of ERF1-V can enhance the salt tolerance of wheat, and heterologous overexpression of GmERF7 can enhance the salt tolerance of tobacco [26,27]. Under alkaline stress conditions, heterologous overexpression of GsERF71 and GsERF6 from wild soybean and VaERF3 from red bean plants can enhance the resistance of Arabidopsis to alkali stress [28][29][30], and overexpression of ZmEREB180 in maize can enhance maize submergence tolerance [31]. Heterologous overexpression of VaERF092 and ERF105 also enhance Arabidopsis cold tolerance [32,33]. Moreover, overexpression of GmERF75 in Arabidopsis can enhance the osmotic stress tolerance of Arabidopsis, and GmERF75 can promote osmotic stress tolerance in transgenic hairy roots [34]. ERF genes can also enhance plant resistance to pathogens. AtERF14 was found to regulate the plant defense response [35]. In Arabidopsis, ERF11 and ERF15 can positively regulate immunity to Pseudomonas syringae [36,37], and in soybean, GmERF13 and GmERF5 can enhance resistance to Phytophthora sojae [38,39]. However, few ERF family genes have been reported to be involved in the response of plants to aluminum stress. In previous studies, the GsERF1 gene was found to be rapidly induced in response to aluminum stress in the wild soybean line BW69 (an Al-resistant Glycine soja line) [3]. However, it is unclear whether GsERF1 is involved in regulating aluminum tolerance in plants. Therefore, the function of GsERF1 was further investigated to elucidate its involvement in the mechanism underlying tolerance to aluminum stress.

Isolation and sequence analysis of the GsERF1 gene
In this study, the full-length cDNA sequence of the GsERF1 gene was cloned from the BW69 line Glycine soja, which is tolerant to Al toxicity. The primers used were designed according to the homologous GsERF1 gene in Glycine max, Glyma09g52900. The GsERF1 gene consists of a complete open reading frame (ORF) of 369 bp, the sequence of which is 99% identical to that of Glyma09g52900 based on the Phytozome genome database, and encodes a protein of 122 amino acids. The predicted GsERF1 protein contains a conserved DNA-binding domain (AP2/ERF domain) of 58 amino acids, which is reported to be the primary functional region. Alignment analysis revealed that the GsERF1 protein was 68 to 97% similar to other proteins encoded by homologous genes and had similar domains (Fig. 1). The analysis of the ERF gene family indicated that GsERF1 is a member of the B-2 subgroup [14]. Many ERF family genes have been reported to have similar functions and play roles in the response to both biotic and abiotic stresses in plants. Phylogenetic analysis showed that GsERF1 and GmERF5 are closely related and belong to the same branch (Fig. 2).

Analysis of GsERF1 expression patterns
Quantitative real-time PCR (qRT-PCR) was performed to assess the transcript levels of GsERF1 in BW69 plants. The qRT-PCR results showed that GsERF1 was constitutively expressed in the roots, stems and leaves. Under aluminum stress, the expression level of GsERF1 in the roots, stems and leaves was significantly increased; this was the case especially in 1-cm-long root tips, in which the expression level increased by 13-fold (Fig. 3c). Over time, GsERF1 was rapidly induced in response to aluminum stress, and the transcripts of GsERF1 reached their maximum level at 4 h, after which the mRNA transcripts of GsERF1 began to decline (Fig. 3a). Under   a GsERF1 expression in the roots exposed to 30 μM AlCl 3 (0.5 mM CaCl 2 , pH 4.5) for 0 to 24 h. b GsERF1 expression in the roots exposed to 0 to 100 μM AlCl 3 (0.5 mM CaCl 2 , pH 4.5) for 6 h. c GsERF1 expression in soybean root apices (R1, 0-1 cm; R2, 0-2 cm), stems (S) and leaves (L) in the absence or presence of Al stress. The values are the means ± SDs (n = 3). The asterisks show significant differences between the control and Al treatments according to Student's t test: *, P < 0.05; **, P < 0.01 treatment with different concentrations of aluminum, GsERF1 transcription increased with increasing AlCl 3 concentration. When the concentration of AlCl 3 was 100 μM, the level of GsERF1 mRNA was 25 times that of the control (Fig. 3b).

GsERF1 is a nuclear protein with transactivation activity
To determine the cellular localization of the GsERF1 protein, its localization was analyzed by expressing a gene encoding a GsERF1-eGFP fusion protein under the control of the CaMV35S promoter in tobacco epidermal cells. An empty vector (pCAMBIA1302-eGFP) was used as a control. As shown in Fig. 4a, the GFP fluorescence was distributed throughout the whole cells, while the GsERF11-eGFP fusion protein fluorescence was visible only in the cell nuclei. These results clearly indicated that GmERF5 is a nuclear localized.
Yeast two-hybrid analysis was used to determine whether GsERF1 could act as a transcriptional activator. The full-length GsERF1 gene was fused to the GAL4 DNA-binding domain and then expressed in the yeast strain Y2H gold to measure transcriptional activation activity; the yeast cells were grown on control medium plates (SD/−Trp) or selective medium plates (SD/− Trp + X-α-gal). Yeast cells containing pGBKT7 plasmids with only the GAL4 DNA-binding domain were used as negative controls. The results showed that only yeast colonies carrying GsERF1 could activate the expression of the reporter gene and cause the colonies to appear blue on the selective medium plate (Fig. 4b).

Overexpression of GsERF1 enhanced plant Al tolerance
To investigate the effect of GsERF1 under aluminum stress, GsERF1 was overexpressed in Arabidopsis to obtain transgenic lines (Fig. S1). Then, three homozygous lines with high expression were selected for phenotypic identification. Under AlCl 3 treatment , the growth of the GsERF1 transgenic plants and WT plants was significantly inhibited, but the root growth of the GsERF1 overexpression (OX) lines was less inhibited than that of the WT plants was. The statistical results also showed that the relative root growth of the GsERF1 transgenic plants was significantly higher than that of WT plants (Fig. 5a&b). Similarly, the fresh weight of GsERF1 transgenic plants was greater than that of the WT plants (Fig.  S2). The proline content of the transgenic plants and wild-type plants increased after aluminum treatment, but the proline content in the GsERF1 transgenic plants was much higher than that in the wild-type plants (Fig. 5c).
To verify the role of GsERF1 in soybean, we generated hairy root lines in which the gene was silenced through soybean hairy root transformation involving RNAiinterference. We detected that GsERF1 mRNA levels were 38 times higher in the overexpression lines than in the WT. However, GsERF1 mRNA levels were 80% lower in the RNAi lines than in the WT lines (Fig. 6a). To visualize activity of GsERF1, the hairy roots of soybean were stained with hematoxylin after treatment in a solution consisting of 50 μM AlCl 3 for 6 hours. The staining showed that the hairy roots of the control displayed stronger staining than did those of the OX-GsERF1 transgenic lines but weaker staining than those of the RNAi-GsERF1 transgenic lines (Fig. 6b). These results indicated that the amount of Al 3+ binding in the hairy roots was the lowest in the OX-GsERF1 lines, while the amount of Al 3+ binding in the hairy roots was the highest in the RNAi-GsERF1 lines. Taken together, the staining results suggested that GsERF1 overexpression can enhance the tolerance of Arabidopsis and soybean to aluminum stress.

Determination of ethylene precursors
To understand the molecular mechanism through which GsERF1 is involved in tolerance to Al stress, several genes from ERF family members were evaluated to investigate their responses to Al stress, and the ACC (an ethylene precursor) content was determined in Arabidopsis. On the basis of GsERF1 expression upregulated in response to ethylene (ET) (Fig. S3), the changes in ACC content in the GsERF1 overexpression (OX) lines and wild-type (WT) Arabidopsis were determined after 10 days of aluminum treatment. The results showed that the ACC content in the GsERF1-overexpressing plants was higher than that in the wild-type plants, while there was little difference in ACC content between the wild-type plants and overexpression plants in the absence of aluminum treatment (Fig. 7a). These results suggested that ET signal transduction may be involved in the aluminum tolerance pathway induced by the GsERF1 gene. In addition, the qRT-PCR results showed that the expression levels of the ACS4, ACS5 and ACS6 genes, which are involved ethylene synthesis, were significantly increased in the GsERF1-overexpressing plants compared with the wildtype plants under AlCl 3 treatment (Fig. 7b, c&d). To   verify whether ethylene is involved in the aluminum tolerance pathway, GsERF1 overexpression lines and wild-type Arabidopsis were treated with ACC and AlCl 3 together. However, the data showed that with the addition of ACC, the root elongation advantage of the transgenic lines no longer occurred in the presence of 100 μM AlCl 3 (Fig. 9) . These results may imply that ethylene can enhance the tolerance of plants to aluminum stress. However, this hypothesis needs to be further verified by subsequent experiments.

Discussion
Aluminum toxicity has a great influence on the roots of plants and directly affects crop yields. Therefore, it is important in theory and in practice to identify new aluminum tolerance genes and determine their functions. The AP2/ERF family is one of the largest families of plant transcription factors, the members of which are involved in many aspects of plant development and responses to multiple environmental stresses [40][41][42][43]. In this study, the GsERF1 gene, which encodes an ERF transcription factor, was isolated on the basis of the aluminum stress resistance gene expression profile of wild soybean line BW69 [3]. Sequence analysis showed that the GsERF1 protein has a highly conserved AP2 domain with characteristics typical of those of domains of the members of the B-3 subgroup of the ERF superfamily ( Fig. 1) [14]. Furthermore, like many other ERF transcription factors, the GsERF1 protein localizes to the nucleus and has selfactivation activity (Fig. 4). Therefore, we speculate that GsERF1 may play a role in plant biological processes.
According to previous studies, members of the AP2/ ERF family have important functions in plants in response to various environmental conditions and at different growth and development stages. This is the case for various members of the B-3 subgroup: OsERF48 can enhance the tolerance of plants to drought and salt stresses [19], AtERF15 can promote the positive regulation of the immune response in Arabidopsis [37], AtERF096 can reduce the water loss rate in Arabidopsis [44], NtERF5 enhances resistance to tobacco mosaic virus, and GsERF16 regulates plant tolerance to bicarbonate in Arabidopsis [28]. ERF genes also play a key role in plant growth and development. For example, MdERF1B can regulate the biosynthesis of anthocyanins and procyanidins in applie [45], and OsERF1 significantly affects the growth and development of transgenic Arabidopsis [16]. However, no ERF family genes have been reported to be involved in the response to Al stress in plants.
In the present study, GsERF1 was rapidly induced in response to aluminum stress, with the greatest amount of transcription occurring under AlCl 3 treatment (Fig. 3).
This result suggested that GsERF1 may play a certain role in the response to aluminum stress.
In recent years, there have been an increasing number of reports in which soybean hairy root transformation was used to verify the function of soybean genes in tolerance to aluminum stress. Among these studies, hematoxylin staining of the GmME1-OE hairy root tips was found to be lighter with decreasing aluminum content. Moreover, overexpression of GmGRPL resulted in obvious aluminum resistance in transgenic lines compared with controls. Under aluminum stress, the aluminum contents in the hairy roots of GmGRPL overexpression lines were lower than those in the controls, and the former presented higher antioxidant capacity. All the above studies are based on the colorimetric properties of hematoxylin, and the degree of staining of hematoxylin in plants roots increased with increasing degree of toxicity of aluminum ions (acidic conditions). Therefore, visual evaluation of stained regions can be used to detect the accumulation of Al in plant roots [46]. In our study, hematoxylin staining was also used to verify whether the GsERF1 gene functions in the tolerance to aluminum under acidic conditions. The results showed that the hairy roots of the GsERF1-overexpressing lines were less stained than were those of the wild type after aluminum treatment, while the GsERF1-RNAi lines were more stained than the wild type (Fig. 6). This indicates that the hairy roots of the GsERF1-overexpressing transformants have a lower aluminum content than those of the wild type do, while the GsERF1-RNAi lines have the highest root aluminum content. Previous studies have shown that root tip elongation is inhibited when plants are subjected to aluminum stress, so root tip elongation is one of the indicators of aluminum tolerance [47]. A similar phenomenon was also revealed by our results: the relative root growth of the GsERF1 transgenic lines with a higher degree of proline accumulation was greater than that of the wild type (Fig. 5). Previous studies have indicated that the proline content increases after plants experience stress, while a large amount of proline can help plants reduce the damage caused by stress. Therefore, on the basis of our experimental results, we speculated that the GsERF1 gene can enhance the tolerance of transgenic Arabidopsis to aluminum stress.
According to previous studies, plant hormones are involved in the response to stress. When plants are under stress, various hormones respond, and different hormones may interact to form a network of mutual exchange to resist external pressure [48]. In our study, we found that when GsERF1-overexpressing plants were subjected to aluminum stress, they produced more ethylene than the wild-type plants did. Specifically, transcription of the ethylene synthesis-related genes ASC4, ASC5 and ASC6 increased (Fig. 7). These results are similar to those found in other studies in which the ASC1, ASC2, and ASC5 genes associated with ethylene synthesis increased under NaHCO 3 stress in GmERF7-overexpressing lines [28]. Protein phosphatase 2A can reduce the toxicity of cadmium by regulating ethylene production in Arabidopsis, and ASC2 and ASC6 were found to be upregulated under cadmium stress [49]. In the present study, the root elongation advantage of transgenic Arabidopsis no longer occurred when ACC was added to the aluminum treatment solution (Fig. 9). We speculated that, within a certain range, ethylene levels may promote root elongation. However, the speculation by which endogenous ethylene levels within a certain range in plants under aluminum stress could have a positive effect on plant growth needs to be further verified. In addition, the content of abscisic acid was also measured to explore its potential role in regulating the response to aluminum stress. The results showed that there was little difference in ABA content between the GsERF1-overexpressing lines and wild-type Arabidopsis plants (Fig. S4). However, some genes involved in the abscisic acid signaling pathway were found to exhibit varied expression levels in the GsERF1-overexpressing lines (Fig. 8). Under aluminum stress, the transcripts of ABI1 and ABI2 in the GsERF1overexpressing lines were significantly lower than those in the wild type (Fig. 8). Previous studies have shown that ABI1 and ABI2 play key roles in abscisic acid signal transduction and act as negative regulators in abscisic acid signal transduction [50]. The transcription of ABI4 and ABI5, which are positive regulators of abscisic acid signal transduction, in the GsERF1-overexpressing lines was higher than that in the WT. [51,52] Furthermore, RD29B was significantly upregulated in the GsERF1-overexpressing lines (Fig. 8). RD29B is mainly involved in drought, salt stress and abscisic acid responses through independent abscisic acid pathways, resulting in higher permeability and stress tolerance of plants [53,54]. These results suggested that the GsERF1 gene may regulate plant tolerance to aluminum stress through the ET pathway and/or the interaction between ethylene and abscisic acid.

Conclusions
In summary, the GsERF1 gene, which encodes a soybean ERF transcription factor, was induced in response to aluminum stress, ethylene and ABA, and its expression was greatest in the roots. GsERF1 overexpression enhanced the tolerance of the transgenic plants to aluminum stress. Molecular analysis indicated that enhanced resistance to aluminum stress might result an increase in the transcription of ACC biosynthesis-related genes and ABAresponsive and/or stress-responsive genes together with proline accumulation. Therefore, the results suggested that GsERF1 may enhance tolerance to aluminum stress mainly through an ethylene-mediated pathway.

Plant material and stress treatments
Seeds from the wild soybean line BW69 were sown in a growth chamber maintained at 28 °C/25 °C under 70% relative humidity and a 14 h light/10 h darkness photoperiod. BW69 is an Al-resistant Glycine soja line, identified by Zeng et al. in our laboratory [3]. The seeds were germinated in vermiculite, and seedlings of uniform growth were selected and cultivated in nutrient solution (pH 5.8) for three days; the solution was renewed daily. After 3 days, the seedlings transplanted into 0.5 mM CaCl 2 (pH 4.5) solution were evaluated for their response to aluminum treatment [55]. For ethylene stress, hydroponically grown seedlings were placed in an airtight plexiglass chamber, and ethylene gas was released after 2 ml of 40% ethephon and 1 g of NaHCO 3 were dissolved in 200 ml of H 2 O [23].
Seeds of Arabidopsis ecotype Columbia (Col-0) plants were germinated and grown in a growth chamber with the following conditions: 22-24 °C, 60% relative humidity, 100 mol photons m − 2 s − 1 light intensity, and a 16 h light/8 h darkness photoperiod. The seeds were sown and grown in potting media from germination to harvest. For analysis of gene expression in Arabidopsis, seeds were sown on 1/2-strength MS agar plates in darkness for 4 days at 4 °C and then placed in a growth chamber. After 10 days, whole plants were collected as samples, and each sample consisted of at least 10 plants.

RNA isolation, cDNA synthesis and quantitative real-time PCR
Total RNA was isolated using TRIzol (Tiangen Biotech, Beijing, China), and cDNA synthesis was performed by the use of a PrimeScript RT Reagent Kit (Takara). QRT-PCR analyses were performed using SYBR Premix ExTaq ™ II Mix (TaKaRa, Shiga, Japan). Actin3 (GenBank accession No. NP_001276160.2) was used as a reference in wild soybean. The Arabidopsis housekeeping gene actin (GenBank accession No. NP_188508.1) was used as a reference in Arabidopsis. The data were analyzed with the 2 -ΔΔCT method [56], and the primers used for qRT PCR are listed in Supplementary Table S1.

GsERF1 gene isolation and sequence analysis
The GsERF1 gene was isolated from wild soybean line BW69. The full sequence of GsERF1 was amplified via PCR in conjunction with the following primer pair: 5′ -GGA TCA CGC CTC AAGTT − 3′ and 5′-CGA ACC CTA AAT CAT CAG − 3′. The PCR products were inserted into the multiple cloning site of a pLB vector (Tiangen Biotech, Beijing, China), and the positive clones were sent for sequencing. Multiple sequence alignment analysis was performed using DNAMAN software. Homology analysis of GsERF1 and the other 44 reference ERF superfamily genes was performed using MEGA 6.0 software through a neighbor-joining method. The amino acid sequences were obtained from GenBank (http:// www. ncbi. nlm. nih. gov/ genba nk/) and Phytozome (http:// phyto zome. jgi. doe. gov/ pz/ portal. html).

Subcellular localization analysis
To analyze the subcellular localization of the GsERF1 protein, full-length GsERF1 was inserted into the NcoI/SpeI sites of a pCAMBIA1302 vector to generate a GsERF1-eGFP construct. The pCAMBIA1302-GsERF1-eGFP fusion construct was subsequently transformed into tobacco epidermal cells. After 2-3 days, the green fluorescence signals in tobacco epidermal cells were observed under a confocal laser-scanning microscope (Olympus FluoView FV1000, Japan). The excitation wavelengths used were 488 nm for eGFP and 580 nm for RFP, and the resolution was 600 dpi [57].

In vitro transcriptional activation assays
For transactivation assays, the full-length GsERF1 gene was inserted into the EcoRI/BamHI sites of a pGBKT7 vector. The pGBKT7-GsERF1 construct was then transformed into yeast strain Y2H gold, and the transformants were grown on SD/−Trp media (Clontech) at 30 °C for 3 days. After selection of the yeast transformants carrying the GsERF1 gene on SD (−Trp) media, they were transferred to SD (−Trp, X-α-Gal) media to evaluate their transcriptional activation. An empty plasmid was used as a negative control.

Arabidopsis transformation and soybean hairy root transformation
Arabidopsis ecotype Col-0 was used for transformation. The full coding region of GsERF1 driven by the CaMV 35S promotor was inserted into the plant expression binary vector pTF101.1, yielding pTF101.1-GsERF1. The construct was subsequently transformed into Agrobacterium tumefaciens strain GV3101, and then the target gene was transferred into Arabidopsis plants by the floral-dip method [58].
Five-day-old seedlings with unfolded cotyledons were used for soybean hairy root production. For the RNAi construct, 233 bp of the GsERF1 coding region was cloned and inserted into a pMU103 vector. The overexpression vector and RNAi interference vector were then transferred into A. rhizogenes strain K599, after which the plants were transformed with the cells by hypocotyl injection [59]. An empty pTF101.1 plant expression binary vector was used as a control.

Hematoxylin staining
The expression of the GsERF1 gene in the hairy root lines was analyzed, and appropriate hairy root lines were selected for subsequent experiments. The hairy roots were treated with 0 or 25 μM AlCl 3 (0.5 mM CaCl 2 , pH 4.5) for 6 h. After AlCl 3 treatment, the hairy roots were washed three times with sterilized water and then stained with hematoxylin. The dyed roots were subsequently washed in sterile water for 30 minutes, after which they were observed and imaged through a Leica S8APO stereomicroscope (Leica, Germany) [46].

Phenotypic analysis of Arabidopsis tolerance to aluminum stress
To analyze the phenotypes of GsERF1-overexpressing (OX) and wild-type (WT) Arabidopsis under aluminum stress, seeds of T 3 GsERF1-overexpressing and WT plants were used. Among them, three transgenic lines with high expression levels were selected. The seed surfaces were sterilized with 10% sodium hypochlorite for 10 minutes and subsequently washed with deionized water. The sterilized seeds were grown on 1/2-strength MS agar plates in darkness for 4 days at 4 °C. Then, the plates were oriented upright and placed in a growth chamber at 22-24 °C, a 60% relative humidity, a 100 μmol photons m-2 s-1 light intensity, and 16 h light/8 h darkness photoperiod. Seedlings with a root length of 1 cm were selected and transferred to 1/2-strength MS agar media (pH 4.5) with different AlCl 3 concentrations. After 10 days, the length from the base of the rosette leaf to the tip of the taproot was measured with a ruler, and images were taken with a Canon EOS 750d camera [60].

Physiological index assays
GsERF1 overexpression and WT lines were treated with or without aluminum for 10 days, and whole plants were selected as samples. The free proline content was measured as described in detail previously [23]. The ethylene precursor (ACC) and abscisic acid contents were determined using an enzyme-linked immunosorbent assay (ELISA) [61].

Statistical analysis
All experiments involving each group were performed at least in triplicate. The data are reported as the means ± SDs. All the data were analyzed via t tests by Graph-Pad Prism 6.01 software to assess significant differences between the means.